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I. INTRODUCTION 

The construction of effective Hamiltonians has proved 
very useful in investigating complicated physical systems, 
since in this way the number of degrees of freedom can 
be reduced. Conventional methods for such a construc- 
tion are based on canonical transformations or projection 
technique. Often the construction is carried out in an 
approximate way, usually by second-order perturbati 
theory. The transformation of the Hubbard modelE. 
into the t-J modecTQ is a well-known example. 

The method of canonical j-transformations was intro- 
duced by Schrieffer and Wolfifl for transferring the Ander- 
son HamiltonianB into the Kondo modela. This approach 
can be sketched as follows. For a given Hamiltonian of 
the form 

n^Ho + eHi (1) 
one searches for a generator S of the transformation 



n' 



He 



eS 



(2) 



under the condition that the new Hamiltonian Ti' does 
not contain linear terms in e. It can be shown that 
Til = ['5', Hq] is the desired conditional equation for the 
generator S. If this is solved for S and inserted into (||), 
one obtains the second-order result of TC'. Since this con- 
struction involves commutation operations with the orig- 
inal Hamiltonian H, the effective Hamiltonian TC' scales 
with the size of the system. The size consistency is impor- 
tant, since otherwise the results for extensive quantities 
like ground-state energy or magnetization would prove 
inconsistent. For a size-consistent calculation of higher- 
order terms, this simple method must be extended, ilihis 
can be done either by a step-by-step transformationulor 
by the method of continuous unitary transformationsEJ. 
However, these approaches are very formal and therefore 
not transparent. l— ■ 

In contrast, the conventional projection methodlla is 
based on a division of the Hilbert space into two sub- 
spaces: the Up subspace in which one is interested and 
the complementary Ug subspace to be projected out with 



the respective projection operators P and Q — 1 — P. 
Using the projection operators P and Q the Schrodinger 
equation is split into two parts, and the states of the 
Uq subspace can be eliminated. Thus one obtains a 
Schrodinger- type equation (H°p — i?) jV^p) = for the 
Up subspace where \ipp) denotes a state of the Up sub- 
space. The effective Hamiltonian is given by 



nf = PHP - PHQ 



1 



QHQ-E 



QHP. 



(3) 



If the two subspaces interact only weakly, the states of 
the subspaces may be separated by a typical energy A. 
In this case one selects the projection operators P and Q 
for a Hamiltonian of the form in such a way that the 
unperturbed Hamiltonian Tip cannot provide transitions 
between the subspaces. By replacing QHQ — E by A, one 
obtains the second order in e of the effective Hamiltonian 
H'^p. For a size-consistent computation of higher orders it 
is not sufficient to expand Eq. (^) for small QHiQ. In ad- 
dition the energy E and the state \ipp) must be expanded 
consistently. The resulting equations for the separate e- 
orders contain eigenvalue corrections of different ordets. 
Therefore these equations must be solved successivelyll3, 
and it is difficult to obtain closed expressions for higher 
orders. 

In contrasti-tOpthc usual projection method the cumu- 
lant approachliJiij preserves size consistency of extensive 
variables and is therefore a suitable tool for the construc- 
tion of effective HaHiiltonians. A previously developed 
cumulant approachEil can only be applied to systems con- 
sisting of two interacting subsystems. This method is 
generalized here for arbitrary systems. 

The paper is organized as follows. In Sec. || the gen- 
eral cumulant approa ch i s presented. As a test for the 
derived method in Sec. HI the effective Hamiltonian up to 
the fourth order in t/U for the half- filled Hubbard model 
is calculated. The conclusions are presented in Sec. IV. 



Finally, a detailed discussion of generalized cumulants is 
given in the Appendix. 
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II. CUMULANT APPROACH 



TV 



The cumulant approachE3~ll3 has estabUshed itself as 
a powerful technique of many-body theory, which makes 
the investigation of static and dynamical ground-state 
properties of weakly and strongly correlated systems jpoai 
sible. It is known from classical statistical mechanical3'E£l 
that size consistency is attained by expressing extensive 
quantities in terms of cumulants, i.e. a cumulant expres- 
sion for an extensive variable scales with the size of the 
system independent of further approximations. In the 
standard diagram technique size consistency is ensured 
in ami approximation by considering linked diagrams 
only.Ej But a diagrammatic approach is usually based on 
Wick's theorenJH3 which is only applicable, if the domi- 
nant part of the Hamiltonian is a single-particle operator. 
Therefore, the diagrammatic description is restricted to 
weakly correlated systems in which the electron-electron 
interaction may be treated pertubatively. 

In the following we propose a cumulant approach for 
the construction of effective Hamiltonians for arbitrary 
systems. The method presented here is based on a per- 
turbational approach for the Hamiltonian Ti. = Hq + Hi . 
The Hilbert space of the unperturbed Hamiltonian Tio 
is split into the low-energy part Up and the high-energy 
part Uq with the respective projection operators P and 
Q = 1 — P. The states of these two subspaces are as- 
sumed to be separated by an energy difference. Now we 
want to construct an effective Hamiltonian for the Up 
subspace. For that purpose the expression for the effec- 
tive Hamiltonian from Ref. Bll 
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(4) 



must be discussed first. Here (3 is the inverse tempera- 
ture. In this case a system with two subsets of degrees of 
freedom a and b was considered, i.e. the Hilbert space of 
the system is a product space of a and b. In Eq. Zb 
and trfc denote the partition function and the trace of the 
subspace b. Eq. (Q) has been derived from the partition 
function of the whole system by separation of the trace. 
Therefore the recently introduced cumulsttJ|t method for 
the construction of effective HamiltoniantO is limited to 
systems on product spaces. 

Now Eq. ^ is generalized by replacing the trace by a 
projection onto the relevant Up subspace 



^-ipinl 



(5) 



where {■ ■ ■ )p denotes P{. . . )P. This generalization is the 
key step in the present construction. In Eq. a part of 
the degrees of freedom is integrated out by the trace in 
the subspace b. In anlogy to this, insignificant degrees of 
freedom are removed in Eq. (|^) by a P projection onto 
the Up subspace. 

In order to transform Eq. (|^) into a cumulant expres- 
sion we introduce generalized cumulants 



■ N 
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.4 = 1 



N 



lnme«'^" 



=0 Vi 

(6) 



For a detailed discussion of generalized cumulants see 
Appendix 0. We define the operator function 



/(A) -Pin (e-'«)pP 



(7) 

(8) 



where Eq. (H) follows from the definition by using 
series expansions. The expression of the effective Hamil- 
tonian can be rewritten 



-T/Cff 



(9) 



For an actual calculation of the effective Hamiltonian as 
formulated here, it is often more practical to start from 
the Laplace transform F{z) of the operator function /(A) 



F{z) 



(nfp 



n 



z-H 



n 



(10) 

(11) 



It can be shown that the second term of Eq. ( |11| ) does 
not contribute to the zeroth and the first order of Hi and 
therefore 
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(13) 



Now we have obtained the sought after perturbation se- 
ries for the effective Hamiltonian and the separate cumu- 
lant expansions corresponding to the perturbation orders 
in Hi ■ Note that the derived cumulant expression is exact 
and applicable at arbitrary temperature since no restric- 
tive assumptions were used. It can be-shown that the 
previously developed cumulant mcthodcj is a particular 
case of our formalism. 

In the following we want to discuss a special situa- 
tion of the method derived above which also provides 
an interesting approximation. If all states of the rele- 
vant Up subspace have the same eigenvalue Ep of the 
unperturbed Hamiltonian Ho, the formalism can be sim- 
plified, because in this case all cumulant expressions of 
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the forms (. . . 'Ho)p and (Wo • ■ ■ )p vanish, and Eq. (|T 
can be rewritten 



Hi 



1 



z-H 



{nfp + ( Hi 



Z — Lq ~ Til y p 



(14) 
(15) 



Lo is the LiouviUe operator with respect to T^o which is 
defined by L^A — [Ho,^] for any operators A. Eq. (15) 
can be shown by transforming the second term of (14) 
into a series of cumulant expressions and using 



{TLi'hC'TLi) p — {TLi{Lq 



UiTHifp. 



(16) 



Note that Eq. ( |l6[ ) is proved, by using the unit operator 
1 = e^^'^ne^'^o, the identity^ e^'^Ae~^'^° = (,KLo+ni) 
and the definition of generafized cumulant expressions 
(^. If we want to restrict ourself to zero temperature, 
the inverse Laplace transform of Eq. (15) is easly per- 
formed and 



(Wo)p+ lim Hi- 



La 



Hi ^/ p 



(17) 



is then obtained for the effective Hamiltonian. In the 
case Til is small relative to Hq, it is straightforward to 
expand the effective Hamiltonian in perturbation series. 
From Eq. (0) it follows that 



,n=0 



Lo 



Hi 



(18) 



As an alternative to perturbation theory, the evaluation 
of the effective Hamiltonian (^5|) can also be done iif-4b£. 
framework of the Mori-Zwanzig projection methodn3'E^ 
as was shown in Ref. ^ In the following we want to 
restrict ourselves to the perturbation result (Ri 



III. APPLICATION TO THE HALF-FILLED 
HUBBARD MODEL 

We now apply the method described above to derive 
an effective Hamiltonian for the strong-coupling limit of 
the half-filled Hubbard model up to fourth order in t/U. 
As it turned out in tlic-past this derivation is highly non- 
trivial; recently Steii£J has used it as attest for Wegner's 
method of continuous transformationsJlil Therefore it is 
reasonable to utilize the cumulant approach presented 
here for re-examination of this expansion. _ _ 

The Hamiltonian of the Hubbard modelcfa on a gen- 
eral lattice in arbitrary spatial dimension is given by 



H = Hq + Hi, 
1 
2 

Hi=tY, A,, ate 



J, a-, 



(19) 
(20) 

(21) 



where c| ^ and Ci^^ are fermion creation and annihilation 
operators for an electron with spin a on site i. t is the 
hopping integral, U denotes the Coulomb repulsion be- 



tween electrons on the same site and ; 



is the 



occupation-number operator for electrons with spin a on 
site i. All informations about the lattice are contained 
in the real hopping matrix Dij. The Hubbard model 
is one of the simplest models which describes Coulomb 
interaction Hq and kinetic energy Hi . 

The expansion of the half-filled Hubbard model around 
the limit of strong coupling, U ^ t, has a long history, 
because this limit makes several simplifications possible. 
Half filling means that the electron number equals the 
number of lattice sites. Before the presentation of the 
Hubbard model, AndersorO showed that such a sys- 
tem is described by an effective spin, Hamiltonian, the so 
called Heisenberg antiferromagnettJ. This Hamiltonian 
is in accordance with the second order {t/U) expansion 
of the half-filled Hubbard model. 

The calculation of highcr_Drders is often based or 
perturbation theory of Katoo. Thus Klein and Seitzt 
obtained the sixth-order spixL interaction for the linear 
chain, and Bulaevski's resultEj of the fourth perturbation 
order for the half-filled Hubbard modeLixi more than one 
dimension was corrected by TakahashiEj. 

Another approach for the derivation of effective Hamil- 
tonians is Jaased on unitary transformations. Harris 
and Langet3 used such a transformation to obtain the 
second-order perturbation. A transformation for calcu- 
latinftijifiher orders was introduced by Chao, Spalek und 
OlesoH'LJ. Beyond second order their results are faulty, 
because in accordance with Eq. (H) the transformed 
Hamiltonian is constructed using a low-order approxi- 
mation of the generator S and so in higher ordecaierms 
appear which mix the different Hubbard band£j'E3. A 
correct transformation algorithm for deriving a pertur- 
bation series of the Hubbard Hamiltonian was proposed 
by MacDonald, Girvin and Yoshiokallj. Thus they could 
calculate the fourth order of tha effective Hamiltonian. 
SteincZl used Wegner's methocO of continuous unitary 
transformations to obtain the fourth-order perturbation 
of the effective Hamiltonian. 

We now apply the method described in Sec. || to con- 
struct an effective Hamiltonian for the half-filled Hub- 
bard modelErH. At half filling all states without doubly 
occupied lattice sites have the lowest eigenvalue of the 
unperturbed Hamiltonian Hq. These states form the low- 
energy subspace with associated projection operator P. 
This subspace is degenerate with respect to Ho- Terms 
with odd powers in the perturbation Hi vanish due to 
the particle-hole symmetry. This is in accordance with a 
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general theorem shown by TakahashiS Therefore from 
Eq. (|l|) foUows 
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(22) 



for the effective Hamiltonian, where the terms in the sum 
correspond to the perturbation orders. The first term 
of Eq. (|l8|) vanishes because the eigenvalue of the un- 
perturbed Hamiltonian Tig of the states without double 
occupied lattice sites is 0. 

For the calculation of the cumulant expressions it is 
profitable to decompose Hi 



(23) 



k=l 



into eigenoperators hk of the Liouville operator Lq with 
eigenvalues A^. The searched decomposition 
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can be derived by introduction of Hubbard operators 

(28) 
(29) 



Note that describes transitions from empty to singly 
occupied sites, whereas describes transitions from 
singly to doubly occupied sites. 

The second perturbation order is given by the first 
term in the sum of Eq. (^). Using the identity 
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which can be proved by series expansion^ we obtain 



(30) 
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One calculates the cumulant expression in accordance 
with the definition (^ by using serie s ex pansions (for 
more details see Appendix |^ and Eq. (|A^)). Due to the 



properties of the hk, expressions of the form {hk)p van- 
ish. Only the term with fc = 3 contributes so that the 
sum can be evaluated 
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-— {niHi)p. 



(32) 



Since Hi describes a hopping process, Eq. ( p^ ) has a clear 
interpretation. Due to the Pauli principle, contributions 
only result from different spins on the sites involved in 
the process. Time reversed initial states can be treated 
commonly by a spin summation. Therefore only two pro- 
cesses contribute to the second-order perturbation and 
one can write symbolically 



HeS 



2nd ord. 
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(33) 



for Eq. (|3^). The diagrams in Eq. ( p^ ) are to be in- 
terpreted as hopping processes where the initial state is 
located above and the final state of the process is below. 
The arrows with double lines lt or JJ. indicate spins which 
have been fiipped. 

Using creation and annihilation operators Eq. ( p3| ) can 
be transformed into an operator equation. Thereby the 
lattice site summations, the spin summation and the 
pre-factors of the hopping processes in accordance with 
Eq. (^ have to be considered. The sign of the diagrams 
results from the process execution and the fermionic anti- 
commutation relations. From Eq. (p^) follows directly 
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Using spin operators 



Eq. 



HcB 



|) can be written as follows 
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(34) 



(35) 



(36) 



where cr^.o-' denotes the vector of the Pauli spin matrices. 

For the calculation of the fourth perturbation order in 
t/U , the second term in the sum of Eq. (^2|) has to be 
evaluated. The cumulant expression can be simplified 
again by using Eq. (30) 
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4th ord. 
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(37) 



The cum ulant expression has to be calculated by using 
Eq. (A3). Due to the properties of the hk terms of the 
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form {hk)p vanish and therefore Eq. 
formed into 



n) can be trans- 



Weft 
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(38) 



Now the fourth order should be transformed to a spin 
Hamiltonian too. In accordance with Eq. ( pO[ ) all states 
have the eigenvalue all of the unperturbed Hamiltonian 
Ti.0 where a denotes the number of doubly occupied lat- 
tice sites of this state. Since Hi describes a hopping 
process, Eq. (^8|) can be evaluated transparently. For 
this purpose, it is appropriate to classify the contributing 
processes and to treat these classes in separate manner. 

Note that only connected diagrams contribute to the 
effective Hamiltonian. This fact results from the effect 
of cumulant expressions and secures the size-consistency. 
The processes to be considered can be classified accord- 
ing to the number of involved sites. As an example, the 
contribution of the processes with three involved sites to 
the effective Hamiltonian will be calculated in the follow- 
ing. In analogy to this computation, the processes with 
two and four involved sites can be considered too. 
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[b] 
[c] 
[d] 
[e] 
[f] 
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[h] 
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first hopping process does not influence the results due 
to the lattice site summations. Note that the geometri- 
cal position of the sites i, j, k is fixed by the hopping 
matrix Di j. In accordance with the number of doubly 
occupied sites after two hopping processes, the processes 
can be divided into two groups. In the cases [a]-[d] re- 
spectively [e]-[h] one notices no respectively one doubly 
occupied sites. Therefore the processes [a]-[d] contribute 
to the first and [e]-[h] contribute to the second term of 
Eq. (38). This fact influences the sign of the process 
diagrams. 
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FIG. 2. Schematic representation of the diagrams con- 
tributing to process [a] of Fig. |^. 
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FIG. 1. Processes of the fourth-order perturbation with 
three involved lattice sites. The spin labelling is dropped. 
Note that the geometrical position of the sites i, j, k is fixed 
by the hopping matrix -Dij . 

Fig. |l| shows the processes with three involved sites 
which are possible in the case of half filling. The selected 



FIG. 3. Schematic representation of the diagrams con- 
tributing to process [e] of Fig. W 
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FIG. 4. Schematic representation of the diagrams con- 
tributing to process [f] of Fig. |l|. 
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FIG. 5. Schematic representation of the diagrams con- 
tributing to process [g] of Fig. |l|. 

The processes [a]-[d] respectively [e] and [h] are equiv- 
alent to each other in each case since their individual 
steps are analogous. On account of the Pauli principle, 
contributions only result again for different spins on the 
sites i, j. Since time reversed initial states can be treated 
commonly by a spin summation, only two initial states 
must be examined. The diagrams contributing to the 
processes [a], [e], [f] and [g] are represented in the Fig. 
^, 1^, 1^ and 1^. The product of the pre-factors Dij in ac- 
cordance with Eq. (2^) is invariant with respect to index 
exchanging of the form j ^ k. Therefore such exchang- 
ing can also be practised in individual diagrams. The 
diagrams contributing to the process [g] of Fig. |l| have 
a different pre-factor and are marked therefore in the 
Fig. I by *. The pre-factor of these diagrams can be 
transformed into the form of the other diagrams by an 
index exchange i <^ j. 



Expression (p8|) for the contribution HcS 
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processes with three involved sites to the effective Hamil- 
tonian can be written symbolically by using the discussed 
process classification and incidence 
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The diagrams of Eq. (39) are to be interpreted again as 
hopping processes. Eq. (|39|) can be summarized formally 
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and transferred to an operator equation 
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Thereby the lattice sites summations, the spin summa- 
tion and the pre-factors have to be considered. The sign 
of the diagrams results again from the special process 
execution and the fermionic anti-commutation relations. 
Eq. ([4l|) can be summarized 
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using spin operators (|35|). 

If the processes with two and four involved sites are 
also treated in this form, all contributions to the fourth- 
order perturbation are considered and up to fourth order 
one obtains the effective Hamiltonian 
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This result is valid for arbitrary lattices and dimensions 
since all such information is contained in the hopping 
matrix Di j. The derived fourth-order perturbation of 
the effective Hamiltonian agrees with other results wlu' 
have been obtained using canonical transformationstj' 
or projection techniquetil. 

In the case of the linear chain with nearest-neighbor 
hopping, processes with four involved sites do not con- 
tribute to the fourth-order perturbation. Therefore in 
this G^e Eq. (^ ) coincides with the result of Klein and 
SeitzEI. 



IV. CONCLUSIONS 

In this paper we have presented a cumulant approach 
for the construction of effective Hamiltonians. The size 
consistency of the results is always guaranteed by the in- 
troduction of generalized cupajilant expressions. While 
a previous cumulant methoa£j has been limited to sys- 
tems containing two interacting subsystems, the formal- 
ism presented here can be applied to systems with arbi- 
trary Hilbert space for any temperature. 



We have applied the presented cumulant method to 
the half-filled Hubbard model on a general lattice in ar- 
bitrary spatial dimension for which the fourth-order per- 
turbation expansion of the effective Hamiltonian was cal- 
culated. In the past the fourth-order perturbation was 
controversially discussed and, therefore, it can be con- 
sidered as a reference problem of methods for the con- 
struction of effective Hamiltonians. The discussion of 
the strong-coupling expansion for the Hubbard model 
demonstrates the power of the derived cumulant ap- 
proach, it enables a transparent construction of size- 
consistent effective Hamiltonians for systems with arbi- 
trary Hilbert space. The derived fourth-order perturba- 
tion of the Hubbard model agrees with results of other 
methods. 
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APPENDIX A: GENERALIZED CUMULANTS 

Thft-cmnulant expressions defined by Eq. (j^) are to be interpreted as a generalization of cumulant expectation 
valuea23'L^. In the Eq. (|^), the Xi are operators whose powers of Vi are to be considered in the cumulant expansion. 
P denotes a projector onto a subspace of the system. Generalized cumulants are traced back to the definition (||), 
i.e. forming cumulant expressions of operator functions are to be explained by means of power series. For computing 
generalized cumulants, a function of non-commuting operators is to be differentiated in accordance with Eq. (|^). 
Therefore, these cumulant expressions must be calculated with the help of series expansions. 

From Eq. (O) follows by expanding the exponential function 
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where the properties of projection operators were utilized. The first addend of Eq. ( |Al| ) does not contribute since it 
is still to be differentiated. The residual logarithm has the form ln(l -|- x) and therefore, can be expanded as follows: 
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N 



■ N 

n 



ni,...,niv— 
(ni,...,rijv)#(0,...,0) 



?1 S. 



N 



■^N )p 



E 



(m+mi) 



ni!mi! 



ni,...,niv— 
mi , . ..,771^—0 
(ni,...,niv)5^(0,...,0) 
L (mi,...,mjv)5^(0,...,0) 



?.=0 Vi 



— i — r 



■^N )p 



If one differentiates the sums of Eq. ( A2) one obtains: 

2 E <^(^^l,?^l + wi) . . .(5(i/Ar,njv + toat) 



[JL^ ...A^jp-^A^ ...A^jp-- 



ni ,...,niv— 
mi ,. . .,m Af — 
(ni,...,njv)/(0,...,0) 
(mi,...,mN)/(0,.--,0) 



ni!mi! 



(xr...x;^")p(xr 



■ Ajv jp 



From Eq. (A3) one can obtain special cumulant expressions. For example, one finds: 



{Afp = iA)p, 
{ABfp = {AB)p l-{A)p{B)p - l{B)p{A)p. 



(A2) 



(A3) 

(A4) 
(A5) 

(A6) 
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